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Abstract 

Accurate prediction of airframe/ actuation coupling 
is required by the imposing demands of modern flight 
control systems. In particular, for agility enhance- 
ment at high angle of attack and low dynamic pres- 
sure, structural integration characteristics such as 
hinge moments, effective actuator stiffness, and air- 
frame/control surface damping can have a significant 
effect on stability predictions. Actuator responses are 
customarily represented with low-order transfer func- 
tions matched to actuator test data, and control sur- 
face stiffness is often modeled as a linear spring. The 
inclusion of the physical properties of actuation and 
its installation on the airframe is therefore addressed 
in this paper using detailed actuator models which con- 
sider the physical, electrical, and mechanical elements 
of actuation. The aeroservoelastic analysis procedure 
is described in which the actuators are modeled as de- 
tailed high-order transfer functions and as approximate 
low-order transfer functions. The impacts of unsteady 
aerodynamic modeling on aeroservoelastic stability are 
also investigated in this paper by varying the order of 
approximation, or number of aerodynamic lag states, in 
the analysis. Test data from a thrust- vectoring config- 
uration of an F/A-18 aircraft are compared to predic- 
tions to determine the effects on accuracy as a function 
of modeling complexity. 

Nomenclature 

Acronyms 

ASE aeroservoelastic 

GVT ground vibration test 

HARV High Alpha Research Vehicle 

* Flight Dynamics Engineer. 
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HOM 

high-order actuator model 

LOM 

low-order actuator model 

RFCS 

research flight control system 

TVCS 

Symbols 

thrust vectoring control system 

A, B, B', C, D, D' 

aeroservoelastic matrices 

AF\B 

antisymmetric fuselage first 
bending 

AWFA 

antisymmetric wing fore- aft 

AWIB 

antisymmetric wing first 
bending 

AWIT 

antisymmetric wing first torsion 

A p 

main ram piston area 

a 

servovalve port width 

b ( 

output structural damping 

b s 

attachment structural damping 

bo, fco, rn o 

aerodynamic damping, aero- 
dynamic stiffness, and 
control surface mass 

b i, ki, mi 

damping, modal stiffness, and 
mass 

C c 

+ xk 

Cf 

general flow coefficient 

C p 

actuator fluid flow pressure 
constant 

c v 

main valve flow coefficient, = 

F e 

externally applied force 

Gi 

ram position feedback gain to 
servovalve 

I 

identity matrix 



hr 

ki 

k s 

m p 

TTi s 

N 

n 

Pi 

Pr 

Ps 

Pi 

P2 

qc 

SFIB 

SWIB 

SWIT 

s 

St 

V T 

X 

x, y 

Xfa 

Xi 

x t 

X V 
X 5 

x v 

Xq 

X\ 

0 aero 
Aae.ro 


combined output/attachment 
structural stiffness, 

_ 1 , 1 
~V e + Y s 

output structural stiffness 
attachment structural stiffness 
mass of main ram piston 

mass of main ram body 
assembly 

bulk modulus 

number of hydraulic systems 
operating 

pressure difference across 
piston, = Pi — P 2 

return pressure 
supply pressure 

ram pressure on attachment 
side 

ram pressure on output side 
impact pressure boundary 

symmetric fuselage first 
bending 

symmetric wing first bending 
symmetric wing first torsion 
Laplace variable 
servoactuator model filter states 

total fluid volume in one 
actuator hydraulic system 

nonzero value 

state and output vector 

actuator ram moment arm 

actuator demand from control 
system 

actuator ram displacement 

displacement of piston relative 
to main ram body, — xi — x s 

actuator body displacement 

displacement of main control 
valve 

displacement of control surface 
mass 

displacement of structural modal 
mass 

zero-order aerodynamic model 
fourth-order aerodynamic model 


6 actual control surface rota- 

= -g 

6 control surface rates 

6d demanded control surface 

rotation, = 

r] a displacement and rate, all actua- 

tor states 

r; e displacement and rate, flexible 

body states 

Tjr displacement and rate, rigid 

body states 

Ti servoactuator model filter time 

constants 

cj actuator stiffness frequency 

Introduction 

Agility and performance requirements, as well as 
a quest for a better understanding of aerodynamics 
by computational fluid dynamicists, are stimulating 
the aerospace community to earnestly investigate high- 
angle-of- attack flight. Modern aircraft are being built 
or reconfigured for enhanced controllability and ma- 
neuverability in stalled and poststalled flight regimes. 
The demanding control tasks at these flight conditions 
introduce some elements of uncertainty in the modeling 
of the aircraft dynamics. 

Aeroservoelastic (ASE) dynamics include the cou- 
pling of structural, control, sensing, aero, and actuator 
dynamics. Structural dynamics can be accurately, al- 
though laboriously, modeled and validated with ground 
test data. Control dynamics are well understood since 
they are designed and extensively verified before fly- 
ing. The primary errors in control feedback sensing 
come from airdata measurements. Therefore calibra- 
tion is necessary to quantify the airdata measurement 
effects on the integrity of the system, and corrections 
are then made so that the feedbacks also can be reliably 
modeled. That leaves aerodynamics and actuation as 
the remaining components of the modeling process. 

Standard techniques for computing linear unsteady 
aerodynamic forces assume potential flow conditions 
independent of angle of attack. Unsteady measured 
data are not available and more precise modeling with 
computational fluid dynamics codes is not feasible 
for prediction of aeroelastic phenomena at a multi- 
tude of flight conditions. Vectoring the engine thrust 
is a promising technology to compensate for a loss 
of aerodynamic control power; however, it is accom- 
plished at the expense of becoming another source of 
model uncertainty by introduction of another coupling 
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mechanism. Plume states are difficult to determine for 
a variety of flight conditions, power levels, and exhaust 
conditions such as turning vane displacements, nozzle 
area, and unsteady flow conditions. 

At high angle of attack and low dynamic pressure, 
conventional aerodynamic control surfaces lose their ef- 
fectiveness and are often driven to their limits. Cou- 
pling between the airframe and actuation may become 
the paramount concern rather than the complex aero- 
dynamics. In particular, structural integration char- 
acteristics such as hinge moments, effective actuator 
stiffness, and airframe/ control surface damping become 
important. Actuators are customarily represented by 
low-order transfer functions matched to actuator test 
data, and control surface stiffness is often modeled 
as a linear spring. 1 More realistic prediction of air- 
frame/actuation coupling is motivated by the more im- 
posing demands for agility enhancement. 

The important ASE modeling uncertainties are, 
therefore, predominantly a function of aerodynam- 
ics and actuation. This paper considers the impact 
of these modeling uncertainties on predicting ASE 
dynamics, with emphasis on the high- angle-of- attack 
flight regime of an F/A-18 aircraft modified with en- 
gine exhaust turning vanes for thrust vectoring. As 
such, detailed actuator models are employed, which 
consider the physical, electrical, and mechanical ele- 
ments of actuation. These results are compared with 
ASE analyses in which the actuators are modeled with 
low-order transfer functions, and corresponding con- 
trol surface dynamic stiffness is simulated as a lin- 
ear spring. Ground test data are used to investi- 
gate the differences in accuracy. Since sustained high- 
angle-of-attack flight tends to be limited to low-speed 
conditions, because of pilot tolerance and for struc- 
tural integrity reasons, aerodynamic forces are gener- 
ally low. Hence, if the modal dynamics and the ac- 
tuation are adequately modeled and validated from 
ground-testing, it is conceivable that aeroservoelastic- 
ity is reduced to essentially a servoelastic phenomenon 
for this flight regime (unsteady plume aerodynamics 
are not addressed) . This hypothesis is investigated by 
varying the order of approximation, or number of aero- 
dynamic lag states, in computing the unsteady aero- 
dynamics. Flight test data are used to determine the 
effects on accuracy as a function of model complexity. 

Aeroservoelastic stability analyses have been suc- 
cessfully derived for the X-29A Advanced Technol- 
ogy Demonstrator 2,3 and the F-15 Short Takeoff and 
Landing (STOL) Maneuver Technology Demonstrator 4 
at predominantly low angle of attack, but no 
known published works have addressed the impor- 
tant ASE coupling mechanisms for thrust-vectored 
configurations and/or high-angle-of-attack flight. This 
paper describes the ASE modeling procedures and 


compares simpler low-order models with high order 
models of actuators and aerodynamics of an F/A-18 
aircraft modified and flight-tested for high-angle-of- 
attack research. Analysis results are compared with 
test data and conclusions are drawn as to the accuracy 
of the various models for high-angle-of-attack flight. 

Aircraft Description 

The F/A-18 High Alpha Research Vehicle (HARV) 
thrust vectoring control system (TVCS) aircraft 
(Fig. 1(a)) is a modification of a preproduction twin 
jet engine fighter F-18 aircraft. It was built by the 
McDonnell Douglas Corp. (St. Louis, Missouri) and 
the Northrop Corp. (Newbury Park, California). The 
aircraft has fixed wings without stores, wingtip folds 
for aircraft carrier storage, wingtip instrumentation 
pods augmented with airdata sensing devices in place 
of missile launcher rails, outboard ailerons, inboard 
trailing-edge flaps, and inboard and outboard leading- 
edge flaps. All actuators are redundant linear hydrome- 
chanical servomechanisms except the leading-edge flap, 
which is the only rotary mechanical system. The twin 
vertical tails are angled outward 20° from the verti- 
cal, and horizontal stabilators are attached to the aft 
fuselage. Leading-edge extensions are attached from 
the wing root leading edge along the fuselage to the 
canopy. 

An F/A-18 aircraft was modified at the NASA Dry- 
den Flight Research Facility to perform flight research 
at high angles of attack using thrust vectoring and in- 
corporating control law concepts for agility and per- 
formance enhancement, as well as providing data for 
correlation with computational fluid dynamics solu- 
tions. This vehicle is referred to as the HARV TVCS 
throughout this report. Modifications include the ad- 
dition of vanes in each engine exhaust for thrust vector- 
ing and corresponding ballast in the forward fuselage 
to maintain center-of- gravity location; a research flight 
control system (RFCS) containing the high-angle-of- 
attack flight research control laws; an inertial naviga- 
tion system for angle of attack and angle of sideslip rate 
feedbacks; wingtip angle-of- attack vanes and pressure 
probes; and additional instrumentation for loads, vane 
temperatures, and structural dynamics. Hence, the air- 
craft has significantly different systems and structural 
characteristics than the basic F/A-18. Aeroservoelastic 
characteristics were expected to change considerably 
from the basic F/A-18 because of the reconfiguration. 

The TVCS structural modifications include a vec- 
toring system with three engine exhaust vanes radially 
displaced about each engine axisymmetric nozzle. Dur- 
ing activation a vast majority of the time the thrust is 
directed with two of the three vanes on each engine, 
while the idle vane is removed from the exhaust plume 
because of temperature considerations. The horizontal 


3 



stabilator also was modified by trimming 2 percent 
of the trailing-edge inboard area to provide clearance 
for the aft end of the outboard TVCS vane fairing 
assembly. 

Total weight added by the aft vane assembly is 2200 
lb, and corresponding ballast in the forward fuselage 
just in front of the cockpit amounts to 700 lb. An 
additional 1300 lb were added with the installation 
of the spin chute and subsystem modifications (emer- 
gency batteries, modified radome, steel reinforcement, 
wiring, etc.). The total weight increase from the ba- 
sic unmodified aircraft is approximately 4200 lb. The 
weight increase resulted in an aft cent er-of- gravity shift 
of 4.4 in. for the empty aircraft. The internal fuel ca- 
pacity is 10,620 lb and full takeoff gross weight is ap- 
proximately 40,000 lb. 

The F-18 basic control system was modified for the 
TVCS to permit communication with a RFCS. 5 To 
control the vanes, as well as the aerodynamic surfaces, 
the RFCS is designed to admit feedbacks computed 
from an inertial navigation system. The research flight 
envelope (Fig. 1(b)) is restricted to a region bounded 
by 15,000 to 35,000 ft altitude, and up to a Mach num- 
ber of 0.7. Beyond an impact pressure of 150 lb/ft 2 , 
lateral-directional control is achieved by the basic con- 
trol laws augmented with an additional yaw vectored 
thrust command in the rudder command path. The 
logic to convert a vectored thrust command to six in- 
dividual vane commands is a very complicated func- 
tion of estimated gross thrust, throttle position, nozzle 
position, angle of attack, and altitude, among other 
parameters. 5 

Aeroservoelastic Modeling Procedure 

Finite-Element Structural Model 

The TVCS modifications such as fairing assemblies, 
bell cranks, vanes, and hydraulic actuators were incor- 
porated in the basic F/A-18 finite-element model to 
produce the analysis model shown in Fig. 2(a). The 
STARS 6 computer program was used to perform the 
free vibration analysis for the aircraft configured with 
full and partial internal fuel with gear up. Frequen- 
cies, generalized masses, and mode shapes were then 
compared with ground vibration test (GVT) results. 
Table 1 shows good frequency correlation of the sym- 
metric (Table 1(a)) and antisymmetric (Table 1(b)) 
modes between the GVT and predicted results. Ac- 
curacy was essential since these predictions were re- 
quired in the aeroelastic and ASE analyses. The unique 
modes of interest for this thrust- vectoring configura- 
tion are the vane rotation modes about their respective 
hinge assemblies. For example, the top vane first rota- 
tion mode had a predicted value of 22.10 Hz, which is 
slightly higher than the measured value of 20.10 Hz. 


Hence, test responses exhibited more motion of the 
vane/actuator fairing attachment point than predicted 
for the same input force. 

A vectored thrust mode shape is generated for each 
exhaust plume by displacing each of the three vanes ac- 
cording to a specified vectored thrust angle, assuming 
a constant and equal loading on each vane. The load 
used was derived from the maximum afterburner thrust 
estimate, which would be the most common configura- 
tion for high angle of attack. The effects of all six vanes 
are incorporated in this simplified manner to produce 
a vectored thrust angle, and therefore a control surface 
mode shape for vectored thrust. 

Aerodynamic Model 

Subsonic linear unsteady aerodynamics were calcu- 
lated using the doublet lattice method. 7 Figure 2(b) 
depicts the aerodynamic paneling used in computing 
the unsteady aerodynamic forces for purely oscilla- 
tory motion of each mode over a range of specified 
reduced frequency values. These generalized aerody- 
namic forces are tabulated complex values which are 
valid only for undamped oscillatory motion at low an- 
gle of attack. Linear theory may not be applicable 
at moderate-to-high angle of attack. However, ac- 
curate modeling of hi gh-angle-of- attack aerodynamics, 
whether steady or unsteady, is computationally pro- 
hibitive even with a rigid airplane. In fact, precise 
modeling of vortex dynamics with viscosity and tur- 
bulence for an entire aircraft is still a topic of intense 
research. 

Justification for the application of forced oscillation 
linear theory is based on the following assumptions 
from Refs. 8 and 9: 

• The structural vibrations about a fixed mean an- 
gle of attack are relatively small, so the rela- 
tion between system displacements and motion- 
dependent structural and aerodynamic forces is 
linear. 

• Small-amplitude oscillatory motion has little influ- 
ence on the unsteady pressures resulting from flow 
separation. 

Also, for full separated flow above 30° angle of attack, 
it is assumed that: 

• Aerodynamic forces caused by modal vibration do 
not change appreciably as a function of angle of 
attack, since the pressures are essentially constant 
chordwise for a fixed mean angle of attack. 

• At low dynamic pressure where sustained high- 
angle-of-attack flight is commonly performed, 
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unsteady airloads caused by modal vibration are 

small and can be reasonably estimated with po- 
tential theory. 

The generalized aerodynamic forces are extended 
to the Laplace domain using approximate analytic 
continuation in order to represent aerodynamic damp- 
ing. Reference 10 outlines some methods of approxi- 
mating the generalized force coefficients with matrix 
formulations of rational functions in s, the Laplace 
variable. These are referred to as Pad6 approximates 
in the reduced frequency domain, and are estimates of 
the time lags associated with unsteadiness in the flow. 
Subsequently, a linear time invariant state space real- 
ization can be formulated for each flight condition. 

In this study, the method of least squares detailed in 
Ref. 3 was used to approximate a Pad6 partial fraction 
approximation to the unsteady aerodynamic forces as 
a function of reduced frequency. The number of partial 
fractions used to fit the doublet lattice data is the order 
of the fit. Pad6 fits to the doublet lattice aerodynamics 
are designated fourth order for the high-order aerody- 
namic model. The low-order aerodynamic model ap- 
plies only zero-order static corrections to the steady 
aerodynamics. In order to maintain consistency in the 
calculations, these fits were not optimized using pro- 
cedures as described in Ref. 10. Low-order models are 
preferred from a computational point of view because 
they can be computed much more efficiently and re- 
duce the cost of later ASE analyses. Several stability 
analyses need to be performed at each flight condition, 
and the computational cost increases dramatically with 
the order of the system. Since the actuator models (to 
be discussed) add a considerable number of states, the 
lowest order aerodynamic model with adequate accu- 
racy is desired. 

Since the low-order model simply represents the un- 
steady aerodynamics as a residual static correction to 
the steady aerodynamic data, no aerodynamic states 
are added and the unsteady dynamics are nullified. Al- 
ternatively, the high-order model adds four states for 
each mode used in the analysis, and it is a precise lin- 
ear estimate of the unsteady aerodynamics over a broad 
range of reduced frequency (Ref. 10). Steady aerody- 
namic data from a linearized version of the full aircraft 
simulation database are incorporated in the generalized 
aerodynamic stiffness and damping matrices. 3 This 
incorporation thereby imposes the constraint that the 
aerodynamic forces near zero-reduced frequency match 
the wind-tunnel-based measured data. 

Because of the complexity of the scheduling of indi- 
vidual vane displacements to achieve a command vec- 
tored thrust direction (noted previously), the vanes are 
combined to derive a vectored thrust control mode with 
respect to a specified thrust level. Vectored thrust force 


effectiveness is reduced to the steady component of the 
plume dynamics acting on the vane configuration used 
to derive the mode shape discussed previously. The 
orientation of each vane determines the force transmis- 
sion to the aircraft through the attachment structure. 
These force effectiveness data were obtained from tests 
as described in Ref. 11. Unsteady plume effects on the 
vanes are not modeled. 

The assumption is that all the vanes have the same 
loading as a function of flight condition, on the av- 
erage. Since the control system attempts to idle one 
vane in each plume to avoid overheating at all times, 
the assumption is considered reasonable because at any 
flight condition each vane will travel between limits of 
idle and a deflection necessary to achieve a vectored 
thrust direction. While one vane in each plume is 
idling, the other two need to compensate to maintain a 
commanded thrust vector angle. There is no reason to 
expect that any one vane is a dominant load- carrying 
member compared to the others during a dynamic ma- 
neuver. 

Actuator Models 

The validity of an ASE model depends not only on 
the dynamic modeling of its components but also on 
the coupling between these components. The airframe, 
aerodynamics, actuation, and control system must be 
accurately modeled over a designated frequency range. 
Actuator responses are customarily represented with 
transfer functions matched to test results, and control 
surface stiffness is commonly modeled as a linear spring 
in the formulation of the airframe dynamics. For high- 
bandwidth flight control systems, as well as agility en- 
hancement at conditions such as high angle of attack, 
it is preferable to improve the actuation modeling by 
considering the physical properties of the system and 
its compliance with the structure. Reference 12 empha- 
sizes, by example, that detailed models based on elec- 
trical and hydromechanical elements are desirable for 
accuracy. In this report the concern is with the sensi- 
tivity of actuation coupling mechanisms with the struc- 
tural dynamics, control dynamics, and aerodynamics. 

The F/A-18 control actuation devices are redun- 
dant hydromechanical servomechanisms powered by 
3000 lb/in 2 hydraulic systems. The actuator functional 
characteristics are shown in Table 2. The aileron actu- 
ation operation is used to demonstrate the comprehen- 
sive modeling procedure in this paper, which is derived 
from Ref. 13. 

For a linear analysis of the actuation system, the 
following assumptions are made: 

• Linear valve flow 

• No valve overlaps, underlaps, or leakage 
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• Parallel-sided valve ports 

• Equal fluid volumes on each side of the main ram 
pistons 

• Equal areas on each side of the main ram pistons 

• No steady load on the ram 

• No friction, backlash, or hysteresis 

• No limits encountered 

Small amplitude characteristics are affected by over- 
lap, hysteresis, friction, and backlash. Limits, port 
shape, and ram symmetry affect large- amplitude mo- 
tion. Therefore the model is presumed to be most ac- 
curate for medium-range amplitudes. 

Linear modeling begins by separating the elements 
of the actuation system into a plant, which contains 
mechanical and hydraulic components, and a feedback 
system which contains servovalve and main ram feed- 
backs. The diagram in Fig. 3 represents the F/A- 
18 aileron actuator installation and control system. 
The perturbation of a free control surface couples with 
the modal dynamics according to the damping b x and 
modal mass mi of a particular structural mode through 
the hinge moment. Aerodynamic forces with compo- 
nents of damping fro and stiffness ko , furthermore, in- 
duce aerodynamic hinge moments. The actuator is at- 
tached to the control surface by a lever arm (denoted 
by damping bg and stiffness kg) and to the airframe 
with some mounting assembly (fr 5 and k s ). The mo- 
tion of the piston relative to the main ram body, x p , is 
fed back to the commanded signal x t . Feedback is also 
employed in the servovalve. 

The servovalve dynamics are assumed to be first or- 
der and are combined with an inner loop position feed- 
back. The outer loop ram position feedback controls 
the ram position according to external demand x x from 
the control system. Dynamic actuator stiffness is a 
sum of the airframe attachment stiffness and the con- 
trol surface lever arm attachment stiffness. Damping 
comes from these sources as well as structural mode, 
aerodynamic, and hydraulic origins. External force F e 
is designated as acting on the main mass from a ground 
impedance test, for instance, to check stiffness. 

Derivation of linear actuation models is valid for mul- 
tiple hydraulic systems operating the same actuator in 
tandem or parallel. The systems are assumed identi- 
cal and the output force produced by the ram is mul- 
tiplied by n (usually 1 or 2), the number of operat- 
ing systems. For this analysis, output and attachment 
structural dampings bg and b s are assumed zero, and 
the masses of the main ram body m s and piston m p 
are neglected. Referring to Fig. 3, in matrix form the 
main ram hydromechanical dynamic equations become 
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for a sign convention assuming that actuator exten- 
sion always gives negative hinge moment and control 
rotation. Therefore, for a servovalve command x v or 
applied test force the ram piston is displaced x p 
relative to the ram body. These equations include the 
effects of coupling with a structural mode of general- 
ized mass mi, and aerodynamic coupling of the sur- 
face inertia through the mo terms. The constant C v is 
a composite coefficient which may be interpreted as a 
combination of servovalve flow pressure gain and some 
leakage across the piston head. 12 


By augmenting this main ram model with the ser- 
vovalve and ram feedbacks shown in Fig. 3, for the 
F/A-18 aileron actuator (s* states correspond to re- 
spective time constants r x , and X denotes a nonzero 
value, which is a function of the system parameters), 
the equations become: 
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Given the extensive assumptions for linearity and ig- 
noring other system components such as a dynamic 
leak valve, damper accumulator, or some other pres- 
sure feedback modulation, the model is not expected 
to match measured results exactly. A model for fric- 
tion would generally add some damping. Also, the fluid 
coefficients C p , C v , and Cj are susceptible to error, and 
the fluid bulk modulus N may differ depending on the 
type of operating system. 12 The experience with the 
F/A-18 actuators was that the servovalve port width 
a (or equivalently, the main valve flow coefficient C v ), 
needed to be adjusted to account for variations in valve 
port shape and size. 

Discounting the states resulting from structural and 
control surface inertias, these high-order actuator mod- 
els (Table 3) for the F/A-18 HARV TVCS become 7th 
order for the aileron and turning vanes, 10th order for 
the stabilator, 6th order for the rudder, 7th order to 


describe the leading-edge flap, and 8th order for the 
trailing-edge flap, thereby adding 45 states for each 
set of symmetric and each set of antisymmetric con- 
trol modes. The differences in number of states be- 
tween actuators is predominantly determined by the 
type of servovalve and the inner loop servovalve feed- 
back mechanism. 


Alternatively, low-order models conventionally used 
for ASE analysis are typically second to fourth order, 
with a corresponding frequency used to define the dy- 
namic stiffness between the airframe and control sur- 
face. They are generated by simply fitting a low-order 
equivalent transfer function to measured data, with no 
regard for actuator details. Compliance with the struc- 
ture is provide d by the stiffness frequency, defined as 
uj = v^T^, which is derived from a measured dy- 
namic stiffness kr and surface mass mo- The trans- 
fer functions and stiffness frequencies used to generate 
the low-order actuator models in the current study are 
listed in Table 4 (the leading-edge flap is not shown 
because of its relative insignificance in control author- 
ity). Comparison plots between the analytical high- 
and low-order models, and measured test results, are 
shown in Fig. 4 for the aileron actuator. This is also 
representative of the other actuators. Beyond 10 Hz 
the test data phase plots show some lag not evident 
in the analytical results, most likely caused by non- 
linearity and the errors in the linear model discussed 
previously. 


State Space Formulation 


Detailed procedures for combining the structural, 
aerodynamic, and sensor dynamics for an ASE anal- 
ysis are given in Refs. 3 and 14. The structural model 
includes rigid body, flexible, and control surface modal 
dynamics. Cross inertias are also included. In essence, 
disregarding aerodynamic lag states (for simplicity in 
demonstration), the state equations take the form 


f Vr \ 
Ve 
Vr 

\ Ve / 


= [A] 


/ 

V 


y 



Vr > 

Ve 

Vr 

Ve / 

Vr \ 

Ve 

Vr 

Ve / 


+ [ B 
+ [ D 


fS\ 

W 

/ 6 \ 


D' ] 


\i) 


for rigid body states 7j r , flexible mode states t] e (includ- 
ing the structural mode states x\ shown previously in 
the actuator model), control surface displacements 8, 
and control surface rates 8. Actuator states are aug- 
mented by adjusting this matrix quadruple as 
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The equation for control surface rate 6, shown previ- 
ously, is augmented to this model for each actuator to 
produce 


x = Ax 4- B<5<* 
y = Cx + D 6d 


where 

x = \VaVrVeVrVe$\ 

includes all actuator states r/ a . This aircraft plant de- 
scription is augmented with control system dynamics 
to construct the entire ASE model. 

Results and Discussion 

In this section some open-loop comparisons are used 
to investigate the effects of actuator and aerodynamic 
modeling. The emphasis here is to determine the sensi- 
tivity of the ASE analysis to actuator and aerodynamic 
model order. Analytical loop gains are determined and 
compared for the different model complexities. Finally, 
flight data results are employed to ascertain the rela- 
tive accuracy of the models in terms of modal stability 
and frequency- dependent stability margins. 


show the greatest difference between the models of all 
the responses of Fig. 5, a comparison of the analyti- 
cal rudder to feedback open-loop responses with GVT 
data are shown in Fig. 6. In the rudder command to 
roll rate response of Fig. 6(a), Fig. 6(b) (rudder com- 
mand to yaw rate), and Fig. 6(c), (rudder command to 
lateral acceleration) the high-order model is matching 
the data better than the low-order model since the low- 
order model gain tends to be too high. In general, the 
degree of complexity incorporated in actuator modeling 
is an important factor in the response characteristics of 
the dynamic model. 

Sensitivity to aerodynamic modeling is shown in the 
open-loop plots of Fig. 7, where a flight condition of 
30,000-ft altitude, 30 lb/ft 2 dynamic pressure, and 40° 
angle of attack is used for computing the aerodynam- 
ics. Zero-order and fourth-order unsteady aerodynam- 
ics are employed in models with low-order actuator 
dynamics (LOM) in Fig. 7(a), and high-order actua- 
tor dynamics (HOM) in Fig. 7(b), using the same four 
transfer functions of Fig. 5 in each case. Notice that 
the comparisons are essentially invariant with respect 
to the order of the aerodynamic fit for each of the four 
pairs of transfer functions. In fact, at the modal fre- 
quencies i, the gain values of the modal peaks in Fig. 7, 
computed with aerodynamics, are very close to the cor- 
responding values in Fig. 5, computed without aerody- 
namics. This observation lends support to the premise 
that the aerodynamics, of any order at low dynamic 
pressure, have an insignificant effect on modal stabil- 
ity. Also, the type of actuator modeling used has little 
influence on the relative accuracy of the unsteady aero- 
dynamic modeling procedure. Apparently the order of 
the aerodynamic model does not significantly change 
the modal response at these fuselage- mounted sensor 
locations. 


Open-Loop Comparisons 

Data in Fig. 5 compare four pairs of open-loop re- 
sponses from surface commands to sensor feedbacks 
between models with high- and low-order actuators, 
thereby representing a GVT by excluding aerodynamic 
forces. The comparisons range from being almost iden- 
tical (pitch thrust vector command to pitch rate) to a 
different order of magnitude (rudder command to yaw 
rate), which indicates a lack of consistency in general. 
Commanded pitch vectored thrust is very insensitive to 
differences in actuator modeling, yet commanded rud- 
der to feedback responses generally demonstrate sig- 
nificant discrepancy between high- and low-order ac- 
tuation. An order of magnitude is equivalent to 20 
dB, yet the common closed-loop stability criteria for 
gain stabilization is to attain an 8-dB margin for all 
structural modes. Therefore a 20-dB open-loop dis- 
crepancy may result in erroneous conclusions from a 
stability analysis. Since the rudder to yaw rate data 


However, as indicated in Fig. 8 where the same com- 
parisons are made with aerodynamics computed at an- 
other condition of 30,000-ft altitude, 170 lb/ft 2 dy- 
namic pressure, and 5° angle of attack, there is more 
distinction between models computed with no aerody- 
namic lag states and those computed with four lags. 
These higher dynamic pressure and low-angle-of-attack 
dissimilarities are seen to be up to 1 or 2 orders of mag- 
nitude (pitch thrust vector command to pitch rate in 
Figs. 8(a) and 8(b)), which is very serious. Therefore, 
unsteady aerodynamic modeling becomes more impor- 
tant with increasing dynamic pressure. Again, the type 
of actuator modeling does not affect the character of 
the comparisons resulting from aerodynamic model- 
ing, as seen by the similar discrepancies between the 
“ Aaero ” and “Oaero” lines in Figs. 8(a) and 8(b). Ap- 
parently, actuator/aerodynamic coupling in the form 
of aerodynamic hinge moments does not vary signifi- 
cantly as a function of actuator model. 
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Loop Gain Comparisons 

Gain stabilization criteria specify maximum accept- 
able loop gains for the model to ensure stability robust- 
ness of the actual aircraft dynamics. Sets of analytical 
single-input-single-output longitudinal, lateral, and di- 
rectional loop gains are shown in Figs. 9(a) through 
9(d) for flight conditions at 5°, 10°, 40° , and 70° an- 
gle of attack, respectively. Dissimilarities between the 
models incorporating low- and high-order actuator dy- 
namics are apparent at all conditions. At the low- 
angle-of-attack and higher dynamic pressure conditions 
of Fig. 9(a) (170 lb/ft 2 dynamic pressure) and to a 
lesser extent Fig. 9(b) (100 lb/ft 2 dynamic pressure), 
the impact of aerodynamic and actuator modeling is 
distinct especially in the directional axis. The large 
miscomparisons in the directional axis stem from dif- 
ferences in dynamics seen in the rudder to yaw rate 
responses of Fig. 5 (actuator models), and Figs. 8(a) 
and 8(b) (aerodynamic models), in addition to other 
lateral- directional feedback responses that were not 
shown. In fact, in the directional axis at 5° angle of 
attack, a change in gain of nearly 20 dB appears to be 
caused independently by actuator modeling and aero- 
dynamic modeling, thereby summing to a total of ap- 
proximately 40 dB. In Fig. 9(c) (40° angle of attack) 
and Fig. 9(d) (70° angle of attack), the responses are 
from conditions at the same low dynamic pressure of 
30 lb/ft 2 . The only source for differences between the 
plots in Fig. 9(c) compared to Fig. 9(d) is in the control 
laws. The only cause of discrepancy within each plot 
is in the actuator modeling. The low-order actuator 
models seem to exhibit more of a tendency toward in- 
stability with either aerodynamic model, as evidenced 
by the highest loop gains. 

In the low dynamic pressure and high- angle-of- attack 
conditions of Figs. 9(c) and 9(d), the loop gains show 
that differences resulting from aerodynamic modeling 
are clearly negligible compared to the clear distinc- 
tions resulting from actuator modeling. The loop gains 
computed from the low-order actuator models gener- 
ally give much more conservative stability predictions, 
which could result in needless system modifications for 
safety of flight. 

Damping Comparisons With Flight Data 
Results 

Modal frequency and damping estimated from flight 
test data are compared with analytical results in 
Fig. 10 for the same flight conditions as Fig. 9. The 
symmetric wing torsion mode near 11 Hz was consis- 
tently predicted to be at a slightly higher frequency 
of 12 Hz. In general, the predicted symmetric modal 
dampings and most antisymmetric modal dampings 
seem to agree favorably with the flight-derived results, 
and the damping values are not very sensitive to varia- 


tions in the models. However, there are some notewor- 
thy exceptions. 

At the 5° angle-of- attack condition (Fig. 10(a)), the 
low-order actuator result with four aerodynamic lag 
terms yields a damping level that is too conservative 
(unstable) near 9 Hz (antisymmetric wing first bend- 
ing), yet the same actuator modeling with no aerody- 
namic lag terms gives a modal damping near 9 Hz, 
which is too high. This extreme variability between 
high- and low-order aerodynamic modeling with low- 
order actuator modeling at low angle of attack and high 
dynamic pressure was also indicated in the directional 
loop gain of Fig. 9(a). 

Another anomaly occurs at the 70° angle-of-attack 
condition in Fig. 10(d), where the low-order actuator 
modeling produces damping levels that are too high 
compared with the flight estimates near 9 Hz (antisym- 
metric wing first bending) and 13 Hz (antisymmetric 
wing first torsion), even though Fig. 9(d) indicated loop 
gains greater than 0 dB at these frequencies. It seems 
that combining low-order actuator models with either 
aerodynamic model is producing more phase-sensitive 
predictions since gain stabilization criteria are not met 
even though damping is more than adequate. Low- 
order actuator models only account for very simplistic 
coupling between the airframe and actuators, so it is 
expected that the predicted phase is inaccurate at the 
modal frequencies. In addition, the predicted closed- 
loop damping levels are not robust using low-order ac- 
tuator models. 

Open-loop responses from actuator commands to 
feedback sensors, with no aerodynamics modeled, 
clearly show that actuator modeling may have seri- 
ous consequences for model accuracy. When unsteady 
aerodynamic accuracy is varied from zero order to 
fourth order, ASE response may change appreciably 
in higher dynamic pressure regimes, yet negligibly in 
high- angle-of-attack flight at low dynamic pressure. 
Hence, unsteady aerodynamic modeling has a negli- 
gible effect on ASE stability at high angle of attack 
and low dynamic pressure, whereas actuator modeling 
can be consequential at any flight condition. Low-order 
actuator models indicate a tendency toward unreliable 
behavior when augmented with various unsteady aero- 
dynamic models. 

Concluding Remarks 

The consequences of aerodynamic and actuator mod- 
eling on aeroservoelasticity, with emphasis on high- 
angle-of- attack flight, were investigated to reveal the 
important aspects for stability predictions. Aeroser- 
voelastic (ASE) analyses were performed in which the 
actuators are modeled as detailed high-order trans- 
fer functions and as approximate low-order transfer 
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functions. The impacts of the unsteady aerodynamic 
modeling were examined by varying the order of ap- 
proximation, or number of aerodynamic lag states, in 
the formulation. Test data from a thrust- vectoring 
configuration of an F/A-18 aircraft are compared to 
predictions to determine the effects on accuracy as a 
function of modeling complexity. 

High-angle-of-attack ASE analyses require detailed 
actuator dynamics modeling. Detailed modeling of ac- 
tuator hydraulic and electromechanical properties is 
important for any ASE analysis regardless of flight con- 
dition. However, the coupling mechanisms between ac- 
tuation and structural dynamics are more pronounced 
with the decreased dynamic pressure characteristic of 
sustained high-angle-of-attack flight. Unsteady aero- 
dynamic forces are relatively insignificant for modal 
stability analyses at these flight conditions. Closed- 
loop stability predictions compared with flight-derived 
estimates show that a low-order actuator model may 
produce unreliable stability estimates at any flight 
condition. 
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Table 1. Free vibration analysis results for the F/A-18 HARV. 


(a) Symmetric modes. 



Light weight (2500 lb fuel) 

Heavy weight (full fuel) 

Primary motion 

Predicted 
frequency, Hz 

GVT 

frequency, Hz 

Predicted 
frequency, Hz 

GVT 

frequency, Hz 

Wing first bending 

5.83 

6.08 

5.75 

6.02 

Fuselage first bending 

7.97 

8.19 

7.45 

7.76 

Wing first torsion 

11.71 

11.78 

11.64 

11.80 

Stabilator first bending 

13.77 

13.78 

13.68 

13.63 

Fin first bending 

15.92 

15.72 

15.92 

15.68 

Fuselage second bending 

16.39 

16.20 

15.34 

15.23 

Wing second bending 

18.05 

— 

17.05 

17.01 

Wing fore-aft 

18.86 

— 

18.19 

— 

Vane rotation 

22.10 

20.10 

22.10 

20.10 

Inboard flap rotation 

23.70 

— 

23.50 

“ 

Wingtip torsion 

27.52 

— 

27.47 



(b) Antisymmetric modes. 



Light weight (2500 lb fuel) 

Heavy weight (full fuel) 

Primary motion 

Predicted 
frequency, Hz 

GVT 

frequency, Hz 

Predicted 
frequency, Hz 

GVT 

frequency, Hz 

Fuselage first bending 

7.40 

7.25 

6.86 

6.64 

Wing first bending 

8.88 

8.48 

8.61 

8.33 

Wing first torsion 

12.03 

12.20 

12.02 

12.13 

Stabilator first bending 

13.69 

13.58 

13.57 

13.45 

Wing fore-aft 

15.36 

15.25 

14.99 

15.09 

Fin first bending 

15.85 

— 

15.56 

15.31 

Fuselage first torsion 

19.50 

— 

18.76 

22.00 

Fuselage second bending 

21.83 

— 

21.01 


Vane rotation 

22.10 

20.10 

22.10 

20.10 

Inboard flap rotation 

23.33 

— 

23.02 

— 

Fore-fuselage torsion 

24.36 

— 

23.98 
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Table 2. F/A-18 HARV actuation system functional characteristics. 



Vanes 

and 

aileron 

Rudder 

Stabilator 

Trailing-edge 

flap 

Output force, lb 
Surface extended 

13,100 

15,740 

29,940 

18,070 

Surface retracted 

12,090 

13,880 

27,250 

14,330 

Surface deflection, deg 

70.0 

60.0 

34.5 

53.0 

Piston area, in 2 
Surface extended 

4.71 

5.64 

10.50 

6.16 

Surface retracted 

4.40 

5.07 

9.72 

4.95 

Stroke, in. 

4.38 

1.43 

7.12 

8.12 

Maximum horn arm, in. 
Aileron only 

4.0 

1.43 

12.0 

9.08 

Vane only 

9.2 

— 

— 

— 

Output velocity, in/sec 

6.70 

1.33 

7.40 

2.76 

Surface rate, deg/sec 

100.0 

56.0 

40.0 

18.0 

Loop gain, sec -1 

48 

37 

30 

18 

Dynamic stiffness, Ib/in. x 10 -5 

2.7 

6.0 

2.0 

0.9 

Maximum free play, deg 

0.126 

0.573 

0.10 

0.573 

Hydraulic flow, gal/min 

8.1 

2.1 

19.8 

4.4 


Control surface 
TVailing-edge 
flaps 


Table 3. High-order actuator model data. 


Transfer function 

fi c: y i r\14 N (j+1264)(,s+805) 

^ 1U > (5+15.4)(s+76.2}(3+469)(5+1033)[a 2 +2(.98)C304)j+304 i ][5 3 +2(,011)(817)5+817 2 ] 


Ailerons 

Stabilators 

Rudders 

Vanes 


(n n v 1016 ^ (a+356) 

a iu ) [s^ + 2(.62)(69.l)a+69.r* i ](5 a + 2(.94)(392)j+392'- i ][^+2(.82)(746)j+746 2 ](5+17326) 

( o w ml4\ (j+284)(*+808)[j 3 +2(.07)(72>5+72 2 ] 

A > {s-|-384)(s4-1098)[s 2 +2(.6}(30)s-(-30 :j ] [ s‘ , 4-2(.8){127)*+127 a j [j a +2(.9){335)s+335 2 ] [j 2 +2(. 02) (718)3+ 718] 

(2.6 X 10 13 ) 4 ^ s _|_224) (5+431 )(i+607)[!s 3 +2(.006)(3335)j+3335 2 ] 

M fi v 1 fl ^ — (a+356) 

U-° x 1U J (3+725)[3 2 + 2(.79)(75)3-l-75 :i 1|3^ + 2(.79)(475)j+475 3 l|3 2 + 2(.82)(781)g+781 2 ■ 


Table 4. Low-order actuator model data. 


Control surface 

Transfer function 

Stiffness frequency, Hz 

Trailing-edge flaps 

35 2 

s 2 + 2(0.71)(35)s + 35 2 

4 

Ailerons 

75 2 

s 2 + 2(0.59) (75)s + 75 2 

16 

Stabilators 

30.8 2 

2.4 

s 2 + 2(0.51) (30.8)s + 30.8 2 

Rudders 

72 2 

5^ + 2(0.69)(72)5 + 72 2 

44 

Vanes 

75 2 

s 2 + 2(0.59)(75)s + 75 2 

8 
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(a) F/A-18 HARV in flight with the TVCS installed. 
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(b) RFCS flight envelope. 

Fig. 1 F/A-18 HARV TVCS and RFCS flight envelope. 
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Fig. 3 Example aileron actuator installation and control system. 






Fig. 4 Frequency response of aileron actuator high- and low-order analytical models compared with test data. 
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Fig. 5 Open-loop frequency response comparisons of surface command to sensor feedback between servoelastic 
models with high-order model (HOM) and low-order model (LOM) actuators. 
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(a) Rudder command to roll rate. 



(b) Rudder command to yaw rate. 

Fig. 6 Open-loop frequency response comparisons between analytical servoelastic models using high-order model 
(HOM) and low-order model (LOM) actuators, and ground vibration test (GVT) results. 
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(a) Low-order actuator model. 


(b) High-order actuator model. 


Fig. 7 Open-loop frequency responses of surface commands to sensor feedbacks. Comparisons are between models 
with zero-order (0 aero, dashed line) and fourth-order (4aero, solid line) aerodynamic models computed for 30,000-ft 
altitude, 30 lb/ ft 2 dynamic pressure, and 40° angle of attack. 
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(a) Low-order actuator model. (b) High-order actuator model. 

Fig 8 Open-loop frequency responses of surface commands to sensor feedbacks. Comparisons are between models 
with zero-order (0 aero, dashed line) and fourth-order (4aero, solid line) aerodynamic models computed for 30,000-ft 
altitude, 170 lb/ft 2 dynamic pressure, and 5° angle of attack. 
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(a) 30, 000- ft altitude, 170 lb/ft 2 dynamic pressure, and 5° angle of attack. 

Fig. 9 F/A-18 RFCS TVCS loop gains from longitudinal (top), lateral (middle), and directional (bottom) loops 
using models with variations in actuator and aerodynamic modeling. 
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(b) 10° angle of attack. 

Fig. 10 Predicted structural mode damping compared to flight data estimation results (solid line), as a function 
of angle of attack, for the primary modes of interest with variations in actuator and aerodynamic modeling. 
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(c) 40° angle of attack. 
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(d) 70° angle of attack. 
Fig. 10 Concluded. 
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